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We show that it is possible to determine the locus of Fisher zeroes in the thermodynamic limit for the Ising 
model on planar ( "fat" ) 0* random graphs and their dual quadrangulations by matching up the real part of the 
high- and low-temperature branches of the expression for the free energy. Similar methods work for the mean-field 
model on generic, "thin" graphs. Series expansions are very easy to obtain for such random graph Ising models. 



1. INTRODUCTION 

One of the more remarkable results to emerge 
from the study of various statistical mechani- 
cal models coupled to two-dimensional quantum 
gravity is a solution of the Ising model in field [Q . 
In discrete form the coupling to gravity takes the 
form of the spin models living on an annealed en- 
semble of triangulations or quadrangulations, or 
their dual planar graphs. The partition function 
for the Ising model on a single graph G" with n 
vertices 

Z,i„gie(G",/3,/i) =^e^S<-)'^''^^+"S-"' (1) 

is promoted to a partition function which incor- 
porates a sum over some class of graphs {G"} by 
the coupling to gravity. 
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The solution to the Ising model in j^] proceeded 
by first forming the grand canonical partition 
function 
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and then noting that this could be expressed as 
the free energy 

Z = - log / V(t)i V(p2 exp ( -Tr ^ 
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of a matrix model, where we have written the 
potential that generates 0"* graphs. In the above 
01^2 are NxN Hermitian matrices, c = exp(— 2/3) 
and the N oo limit is performed in order to 
pick out planar graphs. The graphs of inter- 
est are generated as the Feynman diagrams of 
the "action" in equ. (^, which is constructed so 
as to weight each edge with the correct Boltz- 
mann weights for nearest-neighbour interaction 
Ising spins. Since the edges carry matrix indices 
the graphs in question are "fat" or ribbon graphs. 

The integral of equ. can be evaluated using 
the results of [|j to give (when h = 0) 
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where g is defined by 

2„3 



g{z) = ic^z-^ + z 



2. ZEROES 
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The idea that the zeroes of the partition func- 
tion could determine the phase structure of a spin 
model first appeared in Lee and Yang's work 
They considered how the non-analyticity charac- 
teristic of a phase transition appeared from the 
partition function on finite lattices or graphs, 
which was a polynomial 



m n 
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for a lattice with m edges and n vertices, again 
with c — exp(— 2/3), y — exp(— 2ft,). They (and 
Fisher |^) showed that the behaviour of the ze- 
roes of this polynomial in the complex y or c 
plane, in particular the limiting locus as m, n — > 
oo, determined the phase structure. For temper- 
ature driven transitions, in zero external field for 
simplicity, the thermodynamic limit of the free 
energy on some class of lattices or graphs {G"} 
becomes 

F(G°°, /?) ^-J^ dcp{c)Hc - c{L)), (8) 

where L is some set of curves, or possibly regions, 
in the complex c plane on which the zeroes have 
support and p{c) is the density of the zeroes there. 

The general question of how to extract the lo- 
cus of zeroes analytically has been considered by 
various authors, notably Shrock and collabora- 
tors for Ising and Potts models. It was first 
observed in ^ that such loci could be thought 
of as Stokes lines separating different regions of 
asymptotic behaviour of the partition function in 
the complex temperature or field planes. More 
recently, the case of models with first-order tran- 
sitions has been investigated by Biskup et al. 
who showed that the partition function of a statis- 
tical mechanical model defined in a periodic vol- 
ume V which depends on some complex parame- 
ter such as c or y can be written in terms of com- 
plex functions Fi describing k different phases, 
where the various Fi are the metastable free en- 
ergies per unit volume of the phases, SRF; = F 
characterises the free energy when phase / is sta- 
ble. The zeroes of the partition function are then 
determined from the solutions of the equations 

SRFi = ?StFra < ^Fk, Vfc ^ /, m, 
PViQFi ~ 3F™) = TT mod 27r. (9) 

The equations (^) are thus in perfect agreement 
with the idea that the loci of zeroes should be 
Stokes lines, since the zeroes of Z lie on the com- 
plex phase coexistence curves KF; = 5R-Pm in the 
complex parameter plane. 

The specific Biskup et al. results apply to mod- 
els with first-order transitions, but we are inter- 
ested here in an Ising model with a third-order 
transition, so it might initially seem that these 



results were inapplicable. We are saved by the 
fact that when considered in the complex tem- 
perature plane the transition is continuous only at 
the physical point itself (and possibly some other 
finite set of points). 

3. FAT (AND THIN) ZEROES 

The determination of the locus of Fisher zeroes 
for the Ising model on random graphs in the ther- 
modynamic limit using the ideas of the previoxis 
section turns out to be rather straightforward, as 
we now describe. Since we wish to match ^F 
between the various solution branches to obtain 
the loci of Fisher zeroes and F ^ log(5(c)) for 
the Ising model on planar graphs, the equation 
which determines the loci of zeroes in the ther- 
modynamic limit is 

log|5L(c)| =log|5^^^(c)|, (10) 

where the low-temperature solution gL(c) and 
the various high-temperature solutions gHi (c) are 
given by solving g'{z) = in equ. (^ for z. 

The resulting curve is shown in the c plane in 
Fig. 1 where it can be seen that in addition to 
the physical phase transition at c = 1/4, an un- 
physical transition with the same KPZ [|| expo- 
nents appears at c = —1/4. The interior of the 
curve is the ferromagnetic FM region and the ex- 
terior the paramagnetic PM and unphysical "O" 
phases, separated by cuts on the imaginary axis 
which we have not plotted. 

The points plotted in Fig. 1 are generated from 
a series expansion of Z in equ. (|^), which is ar- 
rived at by reverting the expression for g{z) and 
substituting the resulting z(g) into equ. (^. Ear- 
lier work reported in |^ obtained similar results 
at lower orders. The form of the expression for 
Z means that the contributions from each of the 
terms in equ. (||) are proportional [0 , so the full 
series for Z can be generated from ^ log {z/g). 

The loci of Fisher zeroes are highly non- 
universal, and we also show the zeroes on "thin" , 
generic random 0"^ graphs for comparison in 
Fig. 2. These models can be thought of as the 
^ 1 limit of the matrix models, rather than 
the N oo planar limit. Similar methods to 
those discussed above also serve in this case where 



Figure 1. The Fisher zeroes on fat (j)'^ graphs in 
the complex c plane. 



Figure 2. The Fisher zeroes on thin (p^ graphs in 
the complex c plane. 



one has mean- field behaviour |Tl|. For the Ising 
model on thin graphs F ~ logS", where S is the 
saddle point action for either the low- or high- 
temperature phase. The equivalent of equ. (10) 
is then 



|2(l-c)^| = |(l + c)2(l-2c) 



(11) 



Potts zeroes and chromatic zeroes are also acces- 
sible on the thin graphs. 

In summary, we have seen that an analytic de- 
termination of Fisher zeroes for the Ising model 
on both fat and thin random graphs is possi- 
ble, and that series expansions are easily ob- 
tainable. The general form of the solution also 
holds on (planar) random quadrangulations and 
(j>^ graphs, and in non-zero field, so all of these 
can also be investigated. 
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